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' ABSTRACT 

It is demonstrated that clusters of galaxies are able to keep cosmic rays for a time 

> 

qq ■ exceeding the age of the Universe. This phenomenon reveals itself by the production 
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of the diffuse flux of high energy gamma and neutrino radiation due to the interaction 

ON ■ 

' of the cosmic rays with the intracluster gas. The produced flux is determined by the 

p\' 

cosmological density of baryons, Q^, if a large part of this density is provided by the 
intracluster gas. The signal from relic cosmic rays has to be compared with the flux 

& . 

produced by the late sources, which can be considered as a background in the search 
for cosmic ray production in the past. We calculate this flux considering the normal 
galaxies and AGN in the clusters as the sources of cosmic rays. Another potential 
cosmic ray source is the shock in the gas accreting to a cluster. We found that this 
background is relatively high: the diffuse fluxes produced by relic cosmic rays are of the 
same order of magnitude which can be expected from AGN in the clusters. In all cases 
the predicted diffuse gamma-ray flux is smaller than the observed one and the diffuse 
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neutrino flux can be seen as the small bump at E ~ 10 6 GeF over the atmospheric 
neutrino flux. A bright phase in the galaxy evolution can be a source of the relic 
cosmic rays in clusters, revealing itself by diffuse gamma and neutrino radiations. We 
found that the observation of a signal from the bright phase is better for an individual 
cluster. 

Subject headings: Cosmic Rays: diffusion, confinement - clusters of galaxies - radiation: 
diffuse, gamma, neutrinos 

1. Introduction 

The observations of X-ray radiation from the clusters of galaxies evidence that a large fraction 
of the cosmological baryonic density, pl° s is given by the intracluster gas (White and Frenk 1991, 
White 1993, White and Fabian 1995) . More specifically the cosmological density provided by 
the mass of the intracluster gas, Vlf is close to that derived from nucleosynthesis J7^ s . If 0.$ = 1 
one finds from the references above Vlfh?/ 2 « 0.05 to be compared with the nucleosynthesis value 
0.009 < n% s h 2 < 0.024 (Copi, Schramm and Turner 1995). 

This observation was recently used by Dar and Shaviv (Dar and Shaviv 1995a, Dar and 
Shaviv 1995b) for calculating high energy gamma and neutrino radiations from clusters of galaxies 
(see references therein for the earlier works). In contrast to their calculations, where some ad hoc 
assumptions were made (most notably one about universality of cosmic ray flux in the galaxies 
and clusters), we shall use here a standard picture of production and propagation of cosmic rays 
(CR). This approach results in a lower flux and a different spectral index of radiation at high 
energy as compared with the calculations by Dar and Shaviv. Our conclusion about this low flux 
is valid for all present-day sources of cosmic rays considered in the literature, namely, normal 
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galaxies and AGN in clusters, intracluster and accretion shocks. 

However, as we shall show here, the clusters of galaxies are able to keep the accelerated 
particles for cosmologically long times and thus the production of gamma and neutrino radiations 
at present epochs can be due to particles accelerated in the past. The effective confinement of 
cosmic rays in the clusters was recently recognized by Volk, Aharonian and Breitschwerdt (1996) . 
They considered the cosmic rays in the clusters from the starburst galaxies in the past and TeV 
gamma radiation produced at present. 

A powerful source of CR production can be the bright phase in the evolution of galaxies. 
The idea of bright phase was first put forward by Partridge and Peebles (1967) and was further 
developed in many works (see e.g. (Schwartz, Ostriker and Yahil 1975, Ostriker and Cowie 1981, 
Thompson 1989, Ozernoy and Chrnomordik 1976) ). The bright phase is basically connected with 
the fast evolution of the massive stars (Zeldovich and Novikov 1971, Tinsley 1980). This process 
should result in explosions and shocks in young galaxies and clusters, and thus in acceleration of 
CR. One particular mechanism, convenient for calculations, is SN explosions. This was studied by 
Schwartz, Ostriker and Yahil (1975), Ostriker and Cowie (1981)) and most recently, in connection 
with cluster of galaxies, by Volk et al. (1996). 

Clusters of galaxies are suitable for observations of bright phase in case it occures late, at 
redshifts z ~ 1 — 2 when the clusters were formed. These observations can include a search for 
young galaxies due to blue-enhanced radiation and in particular due to gamma and neutrino 
radiation produced by CR in clusters. 

The plan of our paper is the following. In Section 2 we discuss the propagation of cosmic 
rays and their maximum energy. In Section 3 the confusing problem of the spectrum of produced 
radiation will be considered. The calculations of diffuse fluxes of gamma and neutrino radiation 
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due to normal galaxies are presented in Section 4. In Section 5 we discuss the other sources of 
cosmic rays in a cluster, namely accretion shocks and AGN. Finally, in Section 6 we shall consider 
the bright phase of galaxy evolution relevant to radiation from clusters. 

We shall finish this Introduction by a short summary of the properties of galaxy clusters 
relevant for our calculations. 

A fundamental value which characterizes a cluster is the richness, or number of galaxies, N g , 
in a group or cluster. This value typically varies from 10 to several hundreds. The clusters with 
N g > 100 are defined as rich ones. 

The total mass of a cluster, determined from the virial theorem, is called gravitational mass, 
M grav . The cluster gravitational-mass distribution is given, according to Bahcall and Cen (1993), 
as 

n(> M) = 4 x l - 5 (^)- 1 exp(-^) h 3 Mpc' 3 (1) 

where M is the total (gravitational) mass of a cluster within 1.5/i _1 Mpc and 
M* = (1.8 ± 0.3) x 10 14 /i _1 Mq. Further on we shall refer to the cluster with mass M* as 
the representative cluster. The fraction of the gas mass in a cluster (relative to the gravitational 
mass) is commonly estimated as (0.05 — 0.1)/i -L5 (see (Bohringer 1995, Bahcall 1995) ). In a 
recent work (White and Fabian 1995) this fraction is found for a sample of 13 clusters to be 
between 10 and 22 per cent. 

The distribution (|l]) covers the range from M grav ~ 10 12 /i _1 Mq (small groups) to 
5 x 10 °h M@. The representative mass of a cluster, according to distribution (|l]) is 
M* = 1.8 x lO 14 /i _1 M . The space density of clusters with M > is ~ 1.4 x 10~ 5 Mpc' 3 . 

The size of a cluster, R c i, is defined as the distance at which a galaxy is gravitationally 
confined. Thus, this size can be determined by the velocity dispersion of galaxies a ~ (< v^ al >) x / 2 
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as 

R d = GM grav / < v 2 gal > (2) 

and is typically 1.5/i 1 Mpc, or roughly between 2 — 3 Mpc. 

The hot gas is confined to a cluster due to gravitational potential of the system. The virial 
theorem implies the equilibrium temperature of the gas 

3kT eq = m p < Vg as >= GM grav m p /R c i (3) 

The galaxies in a cluster move in the same gravitational potential and thus from Eq's(||) and (||) 
one has 

r s „~ mW (3*0~4xlo'( Iiji ^ 7; ) 2 K (4) 

If the temperature of intracluster gas were much larger than T eq , the gas expands, if much lower - 
- it condenses. 

Clusters of galaxies are powerful sources of X-ray radiation (Lx ~ 10 42 — 10 44 erg/s). This 
radiation is interpreted as bremsstrahlung radiation of hot intracluster gas with temperature 
T ~ (1 — 10) x 10 7 K, which is indeed close to the equilibrium temperature. 

Both spiral and elliptical galaxies are observed in the clusters. The spiral rich clusters have 
more than 50% of spirals, while the spiral poor clusters - - only 30% (Gorbatskii 1988) (see also 
(Dressier 1980)). However, many observers now think that elliptical galaxies are more abundant 
in rich clusters (Colafrancesco, private communication). 

More than 10 — 20% of clusters have in the center a bright powerful cD galaxy with an 
absolute magnitude My < — 20 m . This galaxy usually provides a large part of radio emission of a 
cluster. Other radiogalaxies are also observed in clusters. Thus, one can assume the presence of 
AGN's in clusters. 
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2. Propagation of cosmic rays and maximum energy 

In the sections 2-4 we shall limit ourselves only to cosmic rays produced by normal galaxies 
in a cluster. This gives a conservative estimate of the produced flux, which can be considered 
as the lower limit, and also gives reference calculations: using them we shall discuss possible 
contributions of other CR sources. 

Our consideration is based on two main observations. 

• The diffuse flux is directly connected with the cosmological density of baryons and not with 
the gas density in a cluster. 

• The spectrum of radiation is determined by the production spectrum of CR in the sources 
and not by the equlibrium spectrum of CR in a cluster. 

As a result one can notice that the problem we study, is defined by a very restricted set of 
parameters, namely, by the cosmological parameter f^, by the cluster richness N g (number of 
galaxies), by the size and age of clusters R c i and t c i, respectively, and by the generation function 
of CR in a typical galaxy, Q 9 ^ l R {E). 

Using CR production in our Galaxy as a guide, we can assume the generation spectrum in a 
galaxy to be Q 9 p al (E) = A(E + E )~^ with E ~ 1 GeV . For E > 1 GeV one then obtains 




where L„ ~ 3 x 10 40 erg/s and 2.1 < 7 9 < 2.4 (Berezinsky et al. 1990, chapter 3). Note, that 
using our Galaxy as a reference source in a cluster results in a somewhat overestimated diffuse flux 
from the clusters. As it will be understood from the following consideration, we do not need the 
detailed picture of CR propagation in a cluster for the calculation of fluxes. It is only important 
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for us to demonstrate that the escape time of CR out of a cluster is larger than the age of the 
cluster, i.e. that the particles injected by galaxies are confined within the cluster during its age. 

In the case of normal galaxies the relevant estimates can be made using an assumption of 
diffusive propagation. However, in case of a powerful outburst, as in the case of the "bright 
phase" , the galactic wind can play a crucial role in CR propagation. 

Observations of radio emission from some clusters evidence that relativistic electrons and 
magnetic field are distributed there on the megaparsec scale (see e.g. Hanisch et al. 1987 ). 
Magnetic field with a typical strength of a few [iG and a principal length scale of 10 — 100 kpc 
was detected in many clusters of galaxies by the Faraday rotation of radiation from radio sources 
(Kim ct al. 1991 and Kronberg 1994 ). 

The stochastic acceleration of electrons in situ on Alfven waves in the Coma cluster can 
explain, according to Schlickeiser et al. (1987) , the observed synchrotron radio emission on a 
cluster scale if the diffusion coefficient of GeV particles is (1 — 4) x 10 29 cm 2 / s if the equipartition 
magnetic field is assumed {B ps 2 fiG). 

Taking the diffusion coefficient Dq = 2 x 10 29 cm 2 /s for particles with energy Eq = 1 GeV and 
using the dependence D oc E 1 / 3 , characteristic for particle scattering off the magnetic irregularities 
with the Kolmogorov spectrum, one can estimate the escape time from a cluster as 

(6) 



12 ( Rd \ 2 /2 x 10 29 cm 2 /s\ / E \ -1 / 3 



6D(E) 

Thus, for energies E < 2 x 10 GeV the escape time r esc is greater than the age of the cluster 
T c i, for which we shall use here and everywhere else the age of the Universe to = 0.82 x 10 10 /i^g yr. 
At energies higher than E c ~ 2 x 10 7 GeV the escape time becomes less than the age of a 
cluster, Td, and production of the secondaries becomes less effective. The determination of this 
critical energy is sensitive to the dependence of the diffusion coefficient on energy which is rather 
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uncertain and depends on the spectrum of the random magnetic field (see e.g. Berezinsky et al. 
(1990) , chapter 9). In particular, the diffusion coefficient may be almost energy independent if the 
turbulent perturbations in the intracluster medium consist of random shocks and discontinuities, 
as is typical in the case of supersonic turbulence. Particles with energies considerably higher than 
2 x 10 7 GeV can be confined in a cluster in this case. 

The extreme case is given by the Bohm diffusion with D ~ crn ■ This is the minimum possible 
diffusion coefficient for propagation along the magnetic field lines. The confinement condition 
Tesc > Td in this case holds up to energy about 10 9 GeV . 

The effective confinement of CR in a cluster was also recognized by Volk et al. (1996). 

Another reason for the diminishing of the efficiency of gamma and neutrino production at 
very high energies can be related to the maximum energy E max reached in the process of CR 
acceleration at the source. For the diffusive shock acceleration by a SN blast wave moving in 
the interstellar medium, this energy is known to be ~ 10 5 — 10 6 GeV for protons. SN shock 
acceleration in predecessor stellar winds of red giants and the Wolf-Rayet stars gives 10 7 GeV 
(Volk and Biermann 1988). Ensemble of shocks in OB associations may provide E m&yi ~ 10 8 GeV 
(Bykov et al. 1990 ). 

Finally, from phenomenological consideration of the origin of ultra high energy CR in our 
Galaxy, it is known (Berezinsky et al. 1990, chapter 4, and Bird et al 1995 ) that particles with 
energies up to ~ 10 9 — 10 10 GeV are most probably produced in the Galaxy, and the generation 
spectrum can have a power-law spectrum up to this energy. 

In the following, we shall use the value E max ~ 1 x 10 8 GeV, either due to inefficiency of 
confinement or acceleration. 
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3. Spectrum of radiation 



This Section is just a note clarifying the difference in the radiation spectrum for the cases 

Tesc > Tel and T esc < T d . 

Let us consider a nuclear-thin object, i.e one with column density much less than the nuclear 
path-length. We assume that high energy protons diffuse from a central source. The generation 
function of the source is Q p oc E~ la and the diffusion coefficient D(E) oc E v . 

The number of secondary particles i = 7, v produced per unit time and unit volume for the 
power-law spectra of CR nucleons can be given in terms of the 'production yields Yi (Berezinsky et 
al. 1990) as 



where a pp = 3.2 x 10~ 26 cm 2 is the normalizing cross section for high energy pp scattering, njj is 
the number density of gas in the cluster and n p (E, r) is the number density of protons with energy 
E at a distance r from the source. This density n p {E, r) can be found from the diffusion equation 
as 



qi{E,r) = YiO- pp n H cn p (E,r) 



(7) 




(8) 



The flux of i-secondaries from the whole source can be found by integrating (0) 



over the volume 




(9) 



For the case T esc > r c / the integration is limited by the upper limit 




(10) 



and one obtains 



Qi(E) = J,YiQ p {E)a pp n H cT cl ~ E~^ 



(11) 
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The physical meaning of this formula is straightforward: each high energy proton emitted by 
a source has the same probability of interaction, equal to o- pp nnCT c i. This result is exact and 
actually valid for any way of CR propagation, e.g. a rectilinear one. 

For the other case T esc < t c i, the integration in eq. (||) is limited by r max = R c i and 

Ql (E) = Yi Q p (E)a pp n H c^r ~ E~^) (12) 

The physical meaning of this formula is also evident: the probability of interaction depends on the 
escape time , R^ l /D(E), and thus on the energy of the particles. 

Therefore, spectra ( |TT| ) and ( |T2] ) are different. This difference often causes confusion. In 
particular, Dar and Shaviv (1995a) took the spectral index of gamma-radiation equal to that of 
the equilibrium spectrum of the protons in our galaxy (i.e. Eq.(|l2|)), while in fact it is equal to 
the index of the generation spectrum of protons (Eq,(|i"l|)), 

To generalize this note, we mention that the absorption time of protons (or electrons), when 
essential, can play the role of age. 

4. Diffuse fluxes of high energy gammas and neutrinos due to normal galaxies 

We consider pp-production of high energy gamma and neutrino radiation. The number of 
gamma's or neutrinos produced by one cluster at time t during the interval dt can be given in 
terms of the production yields Y% (see Berezinsky et al. 1990, chapter 8 and Table 1) as 

Ni(E, t)dt = YiQ9 p al (E)N g ta pp n H cdt, (13) 

where i = 7, u, Q 9 p al is given by Eq.(j5|) and for other notations see Section 3. The values of the 
yields are reported in Table 1. Note that due to the definition of the yields, Y{ and Qi are taken 
at the same energy E. 
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Integrating over time and summing the contribution of all clusters one finds the diffuse flux as 

h{E) = ^ Yi ^^ Q9 ^ E)N ^ n ^ (14) 

where we used n c iMg l as = £,flbPcn with p cr = 1.88 x 10~ 29 h 2 g / cm 3 being the critical density of 
universe. 

In Eq.(|l4|) we neglected the effects of cosmological evolution of the clusters. This is due to the 
fact that we limited ourselves here by normal galaxies as the sources of cosmic rays. Cosmological 
evolution of normal galaxies is weak (X-ray observations) or absent (radio-observations). Besides 
the maximum redshift for the clusters {z max ~ 2 — 3) is small to produce the considerable effects 
for CR. The recent detailed calculations (Colafrancesco and Blasi, in preparation) confirm that 
most of the contribution to gamma rays and neutrinos due to clusters comes from z <J 0.1. The 
exceptional case is given by the burst of CR production (bright phase) which will be considered in 
Section 6. 

From Eq's(||) and (|) one can see that the diffuse flux is determined by the cosmological 
density of baryons fib, the cluster parameters N g and R c i and the galactic quantity Q^ al (E). 

We can now compare the calculated gamma-ray flux (jl4|) with the one extracted from the 
EGRET observations at 10 - 50 GeV (Osborne et al. 1994): 

I° bs (E) = 9.6 x W^E^J 1 cmr 2 s- l sr- 1 GeV- 1 . (15) 

To fit the data we choose 7 9 = 2.11, which is in the allowed range of the generation spectrum 
index. Using £ = 1/2, N g = 100, = 0.116 and R d = 1.5/t -1 Mpc we obtain: 

I 7 (£) = 3.4 x 10- 10 E- 2A1 h 8 o^- cm- 2 s- x sr- x GeV- x . (16) 
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The calculated flux is only ~ 3 x 10~ 4 of that observed experimentally. Notice that one cannot 
diminish arbitrarily the radius of clusters, i? c ;, because the mass inside a smaller radius cannot 
provide Q b h 2 = 0.02. 

This result is not surprising, because the column density CR traversed in the intracluster 
gas, x ~ cmnnHto is of order 0.5 — 1 g/cm 2 , i.e. of the same order as in a galaxy. It means that 
fluxes of gamma-radiation produced by CR in the parent galaxy and in a cluster are of the same 
order of magnitude. It is known that normal galaxies cannot produce the observed diffuse flux at 
energies 100 MeV — 10 GeV. Cluster production of neutrinos with the highest possible energies 
is of greater interest, because the galactic grammage x ga [(E) diminishes as E increases, while 
grammage in a cluster does not depend on energy up to E p ~ 10 8 GeV. 

The diffuse neutrino flux calculated with the help of Eq.(|i~4|) is shown in Fig.l together 
with vertical atmospheric neutrino fluxes. For the prompt atmospheric neutrinos we took the 
calculations of Thunman et al. 1995 (see also references therein). For the calculated flux from 
clusters we used the same parameters as above and neutrino yields from Table 1. One can see 
that for the case of normal galaxies, as CR sources, the predicted neutrino flux is practically 
undetectable. 

The upper limit on the neutrino flux shown in Fig.l is obtained from the condition that the 
diffuse gamma-ray flux observed at 10 — 50 GeV ([15]) is produced by clusters. It results in 

rmax ( rp\ Vp+Ufj. jobs / t?\ 

5. Other sources: AGN's and shocks 

In this section we shall discuss the other sources of CR in the clusters, namely AGN's, shocks 
inside the clusters and the accretion shocks. 
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5.1. Active Galactic Nuclei 

Since about 1% of galaxies have an AGN, we can expect one AGN per cluster with richness 
N g = 100. In particular, roughly 1% of spirals are Seyfert galaxies with luminosity 10 44 — 10 erg/s 
(Weedman 1977). About 5% of elliptical galaxies are radiogalaxies (Schmidt 1978). A very 
powerful source of CR's can be the cD-galaxy often observed in the center of clusters. Recently 
Ensslin et al. (1996) found that radiogalaxies in clusters can heat the intracluster gas due to 
injection of CR. 

There are no well elaborated mechanisms of particle acceleration in AGN's. When particles 
are accelerated diffusively by the accretion shocks in the vicinity of the black hole, they are 
dragged by the flow of gas onto the black hole. A plausible mechanism of acceleration is due to 
unipolar induction in the accretion disc around the black hole (Blandford 1976) or acceleration in 
jets (Biermann and Strittmatter 1987, Quenby and Lieu 1989). It can be expected that up to 10% 
of the total energy release is transferred to accelerated particles. 

However, even assuming one AGN with total luminosity Ltot ~ 10 45 erg/s, we arrive 
at Lqr ~ 10 44 erg/s, to be compared with the total CR luminosity of normal galaxies 
N g L 9 CR ~ 3 x 10 42 erg/s. The gain factor, 30, is still not enough to reconcile the Eq.s (|i~5| ) and 
(|l~6|). In the case of high energy neutrino radiation, this factor is important: the neutrino flux from 
clusters can produce a bump over the atmospheric neutrino flux at E ~ 10 6 GeV (see Fig. 1). 

5.2. Shocks inside clusters 

These shocks can be produced by galactic winds. Even if this phenomenon occurs, it is 
energetically provided by SN induced shocks from the galaxies . CR accelerated in a galaxy obtain 
a total energy which is about the same order of magnitude as one carried away by the galactic 
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wind. Actually the situation is even worse: as was pointed out by Breitschwerdt, McKenzie and 
Volk (1991) and Breitschwerdt, McKenzie and Volk (1993), the intracluster pressure is so large 
that the winds from the galaxies do not develop, except in the case of the starburst phase. 

5.3. Accretion shock 

In this work we assume that a large fraction of cosmologically produced baryons are 
concentrated in the clusters of galaxies. It implies a large rate of accretion of baryonic gas to the 
clusters. This rate averaged over the age of Universe can be roughly estimated as 

< M >~ M gas /t (17) 

where M gas is the mass of the gas in a cluster. 

We imploy in this section the calculations of Bertschinger (1985) for spherically symmetric 
accretion of a mixture of non-dissipative gas (dark matter with a critical density) and a dissipative 
(baryonic) gas. The process of fluctuation growth in both linear and non-linear regimes is 
dominated by the non-dissipative component, while baryonic gas moves in the gravitational 
potential produced by dark matter. However, as demonstrated by Bertschinger (1985) both 
components have the same density distribution p oc r~ 9 / 4 and actually the independent solutions 
for dissipative gas with adiabatic index 7 = 5/3, and non-dissipative gas are very similar. The 
main difference is the existence of a shock in the baryonic gas. The shock originates because of 
the collision of the flow with the dense self-produced core. It propagates outward according to 
R s h ~ t 8 / 9 , where R s h is the shock radius. At present, t = to, for intervals of time r <C to we can 
consider the shock stationary one. 

The basic parameter in this theory is the turnaround radius Rt a - 

The formation of a cluster takes place in the background of the expansion of the Universe. R ta 
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can be defined approximately as the distance where the free-fall velocity Vff = \J1GMjr is equal 
to the Hubble velocity vh = Hqt. As an order of magnitude this condition gives Rf a ~ GM/Hq, 
where M is the total mass inside the turn- around radius Rt a - In our calculations we shall use a 
more accurate expression from (Gunn and Gott 1972, Bertschinger 1985): 

R ta = (SGM^/vr 2 ) 1 / 3 . (18) 

The solutions of Bertschinger are scaling ones which are presented in terms of dimensionless 
quantities. As a dimension parameter for our case we can take the turnaround radius with value 
Rta = 5 Mpc. It gives the gravitational mass inside R c i = 2 Mpc, 

M T V = ^ Pcr <M(A d ) = 3.1 x 1O 14 M , 

where p cr is the critical cosmological density, A = rjR% a is the dimensionless distance and M (A) is 
dimensionless scaling mass given by Bertschinger (1985). In particular, for A c / = R c i/Rt a = 2/5, 
M(\ c i) = 3.7 (in this Section we use h = 0.75). According to Eq.(|l|) this mass corresponds to a 
representative cluster for our consideration. 

The shock radius is found by Bertschinger to be R s h = 0.347i?t a . Normalizing the density 
of baryonic gas by the condition M^ s /M^[ av = 0.1 one can find the gas density at R = R sn as 
Pgas = 2.2 x 10 -29 g/cm 3 . The radial velocity of gas at the position of the shock is, according to 
Bertschinger (1985), v r = l.5R ta /t = 8.4 x 10 7 cm/s. 

The particles are accelerated diffusively on the shock front which, as was mentioned above, 
can be considered as stationary at the present time. The particles are scattered off magnetic 
inhomogeneities upstream and downstream the accreting flux and dragged by it inside the cluster. 
The CR luminosity can be calculated assuming that 10% of kinetic energy of the infalling gas is 
transferred to accelerated particles: 
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L CR = 0.2n Pgas (R sh )R 2 sh vf = 2.3 x 10 44 erg/s. (19) 

The maximum energy of accelerated particles can be estimated in the usual way. Acceleration 
is efficient up to energy at which the acceleration length, D/u becomes comparable with the radius 
of the shock front R s h (see e.g. Jones and Ellison 1991 ). Assuming the Bohm diffusion with 
D ~ r^c one obtains: 

9.3 x 10 17 (H/lfiG)(R sh /3 Mpc)(u r /W 7 cm/s) eV (20) 

The alternative solution for gas accretion to the cluster is given by the pancake model 
(Zeldovich 1970). The shock is also produced in this model. The results of numerical calculations 
for the cluster-size pancake are given by Doroshkevich, Zeldovich and Sunyaev (1976). The 
distance to the shock in this simulation is about 0.6 Mpc for a typical size of a cluster 2 Mpc. The 
gas velocity at the position of the shock is 1.5 x 10 7 cm/s. The CR luminosity can be estimated 
very roughly for the given parameters and geometry as ~ 10 43 erg/s. 

6. Bright Phase 

The first-generation stars were produced due to non-linear growth of the fluctuations in 
the baryonic gas. Their mass spectrum was strongly dominated by the large masses with a 
fundamental scale given by the Jeans mass Mj ~ 10 5 Mq. The fragmentation of these massive 
protostars results in the formation of the initial mass function (IMF), being strongly enhanced 
by massive short-lived first-generation stars. Most of these massive stars finish their evolution 
as SNII. The explosions and stellar winds enrich the interstellar gas by heavy elements and this 
composition is seen in metal deficient population II stars. The duration of the longest phase in the 
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heavy-star evolution, the main sequence stage, can be estimated as the ratio of the thermonuclear 
production of energy, eMc 2 , to the luminosity L, or 

™ S ~ (1.3xlQ38M/M ) (L/L ^ ) ~ X 10 — yr ' 

where the conversion coefficient e 3.3 x 10 -3 is taken for a typical presupernova II massive 
star (M ~ 30M Q ) with 2M Q carbon-oxygen core and 12M helium shell. For these stars 
L/Lzdd ~ 0.1 — 0.01 and thus the typical duration of the bright phase is 10 7 — 10 8 years. Partridge 
and Peebles (1967) in their pioneering work assumed the duration as 3 x 10 7 yr. 

The total energy release during the bright phase is estimated from nuclear conversion. 
Partridge and Peebles (1967) found W~3x 10 61 ergs per galaxy; Schwartz, Ostriker and Yahil 
(1975) estimated it as W ~ 10 61 — 10 62 ergs. These authors considered this energy as being 
released in the form of kinetic energy of SN shells, which results in the shocks and acceleration of 
the particles (see also (Ostriker and Cowie 1981)). An emphasis was given to these processes in 
the clusters of galaxies. In the early works the bright phase was assumed to be at relatively large 
red shifts, e.g. 2 < z < 10 in the work by Schwartz, Ostriker and Yahil (1975). 

In most recent calculations of Volk et al. (1996) the bright phase (starburst phase) was 
placed at relatively low red shift z ~ 2 — 3 and the total energy release was found from the iron 
abundance in the clusters. In the discussion below we shall follow this approach. 

Finally we mention the detailed calculations by White and Frenk (1991) for the galactic 
evolution. These calculations start with the cosmological evolution of fluctuations described in the 
Cold Dark Matter cosmological model and proceed with the stellar evolution up to the present 
age of the Universe. The star formation in this model occurs late: half of the stars are born at 
the epoch with red shift z < 1. No bright phase is found. This result is not a surprise, because 
the authors use the mass function as it is observed at present time, while the main assumption 
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which leads to the bright phase is the strong enhancement of the mass spectrum by supermassive 
short-lived stars. 

How much can we increase the calculated diffuse fluxes taking into account the bright phase 
production? 

The increase can be estimated as the ratio Wq R /Wq R , where Wff R and Wq R are the CR 
energy release in a cluster during the bright phase and due to present-day galaxies, respectively. 
The value Wq' r is given by the product of N$n, the total number of SN explosions in a cluster 
during the bright phase, and Wq R = 2.8 x 10 49 erg, the total energy of CR produced by one SN. 
The former can be found, following Volk et al. (1996), as Mp e /0.1M Q , where Mp e is the mass of 
iron in the intracluster gas and O.1M is the mass of iron produced by one SNII explosion. Since 
the mass fraction of iron in a cluster is about e^ e s=s 7 x 10~ 4 we obtain for our representative 
cluster with the mass of gas Mf 3,8 = 1.8 x lO 13 /i _1 M : 

M gas 

w ™ = €f <mm; w ™ = 44 x 1060 erg - (21) 



The value of Wp R is 



W 9 CR = L p t Ng « 7.7 x 10 59 erg (22) 



and thus the factor of the diffuse flux increase is rather small, about 6. 

Note that W^ R as estimated by Volk et al. (1996) is ~ 15 times larger than our value. The 
main reason for this discrepancy is in the value of W^ R . Volk et al. adopted it from theoretical 
considerations as 3 x 10 50 erg. Our value is found from phenomenology of CR in our Galaxy. 
Namely it is given as Lp/usn, where usn ~ 1/30 yr is the frequency of SN explosions in our 
Galaxy and Lp is the CR luminosity of our galaxy, given by Drury et al. (1989) as 
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L P = cM, 




x{E) 



3 x 10 40 erg/s, 



(23) 



where ujcr{E) is the observed density of CR in particles with kinetic energy E, x(E) is the 
column density of gas (" grammage" ) traversed by particles with energy E in the Galaxy and M g 
is the total mass of the interstellar gas. For a reasonable set of these values our present estimates 
coincide with the luminosity given by Drury et al. (1989). 

As a result, our value of « 2.8 x 10 49 erg is 11 times less than that of Volk et al. (1996). 

We conclude, therefore, that the increase of the diffuse flux due to bright phase as compared 
with normal galaxies is smaller than that due to AGN and accretion shocks. 



We have calculated the diffuse fluxes of high energy gamma and neutrino radiation from 
clusters of galaxies. These radiations are produced in pp-interactions of high energy CR 
with intracluster gas. The basic assumption which we use is that a considerable part of the 
cosmologically produced baryons are concentrated inside the clusters. An interesting observation 
is that, with this assumption, a very restricted set of parameters enters the calculations. They are: 
the generation function of CR in the sources, Qcr{E), the cluster size R c i and the richness N g . 
The exact picture of propagation of CR in a cluster does not affect the produced flux of radiation, 
because the escape time of CR from a cluster is larger than the age of the cluster. For the same 
reason, the spectrum of the produced radiation is very flat. At high energy its exponent equals to 
that of production spectrum of CR in the sources. 

We made explicit calculations for CR produced by normal galaxies. In this case the generation 
function, Qcr(E), for a galaxy can be taken as that of our Galaxy. It gave us a reference case for 



7. 



Discussion and conclusions 
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calculations: for other sources the produced flux is scaled as the CR luminosity of a source Lqr- 
What other possible sources could there be? 

The general restriction is that they must be very powerful to exceed the total CR luminosity 
of normal galaxies in a representative cluster, 3 x 10 42 erg/s. 

The anticipated sources could be AGN's (including the central cD-galaxy in a cluster) and the 
shock in the accretion flow of the gas onto clusters. They can increase Lcr and thus the diffuse 
flux by factor ~ 30. According to our estimates the bright phase can increase Wcr, and thus the 
diffuse flux, only by factor 6 as compared with normal galaxies. 

The normal galaxies in the cluster can provide neither the observed diffuse flux of gamma 
rays nor the detectable neutrino flux at any energy. The accretion shock acceleration and/or 
acceleration of particles in AGN's can result in the production of a detectable bump above the 
atmospheric neutrino flux at energy E ~ 10 6 GeV, while the predicted 0.1 — 100 GeV flux of 
gamma radiation remains one order of magnitude less than observed. 

Our results differ from calculations of Dar and Shaviv (1995a, 1995b) : they are more 
pessimistic for gamma radiation and more optimistic for neutrino radiation. 

Dar and Shaviv assumed that CR density in clusters is the same as in our Galaxy. We 
cannot reproduce their assumption studying all possible sources of CR in the cluster. For the CR 
production by the normal galaxies the energy density of CR in a cluster is 

u d R = ^gjVo = 9 3 x 1Q -4 h 2 e y /cm 3 

to be compared with the galactic value 0.5 eV/cm 3 . Here Lp = 3 x 10 40 erg/s is the CR luminosity 
of our Galaxy and N g = 100 is the richness of a typical cluster. Increasing the CR luminosity of a 
cluster up to 10 44 erg/s, due to accretion shock or AGN, we are still left with more than one order 
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of magnitude deficit. 

The bright phase deserves a special discussion. If SNe are the main sources of CR, the 
bright phase gives the increased flux described as compared to normal galaxies by a factor 
n Siv/ n 57V' w bere n b g N and nlf N are the SN space density due to bright phase and normal galaxies, 
respectively. They can be found as 

nb SN = 7^jrr n(M)M gas dM 
V.UVIq JM min 

and 

POO 

n SN = »SNt / n(M)N g (M)dM 

where n(M) is obtained by differentiating eq. ([!]), M is the gravitational mass and M gas ~ 0.1 M 
is the gas mass in a cluster; N g (M) is the number of galaxies in a cluster with mass M, and i>sn 
is the SN frequency in a normal galaxy. The ratio ng N /n^ 9 N is about 8, i.e. practically the same 
as the previously estimated value assuming all clusters having the characteristic mass M*. 

As we demonstrated the spectrum of high energy radiation coincides with the production 
spectrum of CR. That is why our spectrum can be as flat as E~ 2A , while Dar and Shaviv explicitly 
state that their spectrum is E~ 2 ' 7 . This is the reason that our flux of neutrino radiation is higher 
for the same CR luminosity. 

The bright phase can give a much more pronounced effect for the case of individual cluster. In 
particular for the Perseus cluster, discussed by Volk et al. (1996), the mass of the gas is very large 
M ~ 4 x 10 14 M© for the given richness iV ~ 500, and therefore the total number of SN is very 
large as compared with the number of SN produced by galaxies during their normal evolution. 
Such clusters are the best candidates for observations of both, gamma rays and blue-shifted young 
galaxies within. 

Finally, we shall discuss the possibility of gamma-ray detection from the nearest Virgo cluster 
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(19 Mpc). 

The EGRET observations (Fichtel et al. 1994) give the upper limit for 100 MeV gamma-ray 
flux .F 7 (> lOOMeV) < 0.5 x 10~ 7 cmT 2 s^ 1 , which is actually valid only for the central source M87. 
In the models considered here the flux from the extended region occupied by the whole cluster 
can be considerably larger. Using the richness N=2500 and mass of the gas MrjIx 1O 14 M we 
obtain the flux at E 1 > 100 MeV ~ 10 _8 cm _2 s _1 both for normal galaxies and bright phase. 
The gain factor predicted for a general case is lost for the Virgo cluster due to the huge number of 
galaxies in comparison with the gas mass. Unfortunatelly this flux is undetectable for the existing 
instruments. 
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Table 1: The values of the yields multiplied by 10 3 
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Fig. 1. — The diffuse neutrino flux due to the interactions of the CR produced by Normal 
Galaxies and AGN with the intracluster gas. The upper limit refers to the maximum neutrino flux 
correspondent to the observed diffuse gamma ray flux as produced in clusters. As a comparison 
the atmospheric neutrino flux is plotted with the prompt neutrinos taken into account. 



